Intrinsic Spin-Hall Effect in n-Doped Bulk GaAs 
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We show that the bulk Dresselhauss (fc 3 ) spin-orbit coupling term leads to an intrinsic spin-Hall 
effect in n-doped bulk GaAs, but without the appearance of uniform magnetization. The spin-Hall 
effect in strained and unstrained bulk GaAs has been recently observed experimentally by Kato 
et. al. T]. We show that the experimental result is quantitatively consistent with the intrinsic 
spin-Hall effect due to the Dresselhauss term, when lifetime broadening is taken into account. On 
the other hand, extrinsic contribution to the spin-Hall effect is several orders of magnitude smaller 
than the observed effect. 

PACS numbers: 73.43.-f,72.25.Dc,72.25.Hg,85.75.-d 



Recent theoretical work predicts dissipationless spin 
currents induced by an electric field in semiconductors 
with spin-orbit coupling^, Q- The response equation is 
given by jj = a s eijkEk, where jj is the current of the i-th 
component of the spin along the direction j and & is the 
totally antisymmetric tensor in three dimensions. The re- 
sponse equation was derived by Murakami, Nagaosa and 
Zhang 2] for p-doped semiconductors described by the 
Luttinger model of the spin-3/2 valence band. In an- 
other proposal by Sinova et al. [2j, the spin current is 
induced by an in-plane electric field in the 2-dimensional 
electron gas (2DEG) described by the Rashba model[3j- 
The intrinsic spin-Hall effect predicted by these recent 
theoretical works is fundamentally different from the ex- 
trinsic spin-Hall effect |4U5j due to the Mott type of skew 
scattering by impurities .6| . The intrinsic spin-Hall effect 
arises from the spin-orbit coupling of the host semicon- 
ductor band, and has a finite value in the absence of 
impurities. On the other hand, the extrinsic spin-Hall 
effect arises purely from the spin-orbit coupling to the 
impurity atoms. 

Experimental observation of the spin-Hall effect has 
been recently reported by Kato et al. lj in an electron 
doped bulk sample and by Wunderlich et al. in a two 
dimensional hole gas (2DHG) 7]. The 2DHG experiment 
has been analyzed in a previous paper || where it was 
shown that the vertex correction due to potential impu- 
rity scattering vanishes for that particular system. The 
experimental system is also in the regime where lifetime 
broadening due to impurity scattering is much less than 
the spin splitting, thus strongly suggesting an intrinsic 
mechanism of the spin-Hall effect. In the experiment 
of Ref. 0, spin accumulation due to a spin current is 
observed even in the unstrained GaAs where no appar- 
ent spin splitting is observed. The absence of observed 
spin splitting seems to show the absence of intrinsic spin- 
orbit coupling in unstrained n-doped GaAs. This fact 
prompted the authors of Ref. [l| to interpret the ob- 
served spin-Hall effect in terms of the extrinsic mecha- 
nism due to impurity scattering only. In thispaper we 
show that, under close scrutiny, the results of [jj] are con- 
sistent with an intrinsic mechanism. We first show that 



the unstrained GaAs has a Dresselhauss fc 3 spin splitting 
which escapes detection by the method used in 0, 0. 
We then show that this spin splitting leads to a spin-Hall 
current. This therefore explains the observed spin ac- 
cumulation on the edges of the unstrained GaAs within 
the framework of the intrinsic spin-Hall effect. Further- 
more, the observed magnitude is consistent with the the- 
ory, after lifetime broadening due to impurity scatter- 
ing is taken into account. We also predict that the bulk 
Dresselhauss term produces no net uniform magnetiza- 
tion in the sample, this being generated solely by the 
strain terms. In the case of strained GaAs, we compute 
the self-energy correction in the weak spin-orbit coupling 
limit and find a value for the spin-Hall conductivity close 
(enough) to the measured value. The independence of 
the spin current on the crystalographic directions can 
also be explained by the dominance of the fc 3 term over 
the fc-linear terms induced by the small strain. We also 
perform an order-of magnitude estimate and find out that 
the extrinsic spin-Hall effect is seven orders of magnitude 
lower than the clean limit of the intrinsic spin-Hall effect, 
and several orders of magnitude lower than the observed 
experimental value. 

Let us first examine the extrinsic spin-Hall effect. In 
the extrinsic mechanism^, |5(, there is no spin-orbit cou- 
pling in the band structure, and the spin-Hall effect is 
caused by the scattering of electrons by the spin-orbit 
interaction with impurities. The Hamiltonian is given 
by: 
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where V(r) is the impurity potential. The extrinsic spin- 
Hall effect is basically derived from the atomic Mott scat- 
tering 0|, and the important length scale is governed 
by the Compton wave length A c = h/mc. The extrin- 
sic spin-Hall effect has been computed systematically for 
this Hamiltonian [loj . and the order of the magnitude of 
the effect can be estimated to be: 



O 'extrinsic ~ "T~(A c fc_F) fc_F; 



(2) 
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where kp is the fermi wave vector. For the experimental 
system by Kato et al., with kp = lO 8 ™^ 1 and a con- 
duction band effective mass of m = 0.0665m e , Eq. (J2J 



gives (Textrinsic = 1.2 X 10 4 0" 



1 , almost 4 orders of 



magnitude smaller than the observed spin-Hall conduc- 
tance. On the other hand, the intrinsic spin-Hall effect is 
a genuine solid state effect, governed purely by the fermi 
wave vector kp, and the order of magnitude of the effect 
is given by 



— kp 
n 



(3) 



in the clean limit. Therefore, we see that the ratio of the 
two effects is given by 11] 
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Therefore, in distinguishing between the two effects, it is 
extremely important to keep in mind the smallness of the 
dimensionless parameter X c kp. In the literature of the 
anomalous Hall effect, a so called "enhancement factor" 
is sometimes introduced in a rather ad-hoc basis flpl fl^| . 
However, this "enhancement factor" is microscopically 
based on the spin-orbit coupling within the band struc- 
ture, and would necessarily lead to a spin splitting of the 
bands. Therefore, we can safely conclude that if there 
were no spin splitting due to the intrinsic spin-orbit cou- 
pling within the band, the extrinsic spin-Hall effect is far 
too small to explain the experiment by Kato et at in the 
unstrained GaAs. 

Let us now turn to the intrinsic spin-orbit coupling 
within the conduction band. The Hamiltonian of an 
inversion asymmetric bulk (unstrained) semiconductor 
contains a Dresselhauss fc 3 spin splitting term in the con- 
duction band, which can be written as a momentum de- 
pendent magnetic field: 
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where B x - 

B z = 'yk z (k 2 — k' x ). The coupling constant 7 has been 
determined in a number of independent experiments, and 

o 3 

a value of 7 « 25eVA is widely quoted in the litera- 
ture 0, Q, 0, . We must now reconcile this spin 
splitting with the fact that the measurement carried out 
in Ref. [tj does not see any splitting in the unstrained 
sample. In [jj, a spin packet injected at the Fermi mo- 
mentum is subsequently dragged by an external electric 
field E. Experiments are performed along two crystallo- 
graphic directions £7||[110] and £?||[110]. This creates an 
average nonzero particle momentum (k) ~ § Et which 
in turn creates a non-zero average (over the Fermi sur- 
face) internal magnetic field (B). The spin splitting in 
is obtained as a derivative of the averaged (B) with 
respect to the drag momentum (k). Due to the special 



symmetry of the Dresselhauss spin-orbit coupling, this 
procedure turns out to yield a null result, even if 7 is 
finite. Take, for example £7| | [110] , then the momentum 
of a particle injected near the Fermi momentum k F is: 

k = k F + (k), (k) =-— E\\[U0], (k x ) = (k y ). (6) 
m 

To first order in (k), the components (say x) of (B) av- 
eraged over the Fermi surface is: 



(B x ) = (k x ) J ^ [(fcf f - (k F ) 2 - 2k F k F 



(7) 



Since the spin-orbit coupling term is much smaller than 
the kinetic term, the Fermi surface is, to first order in 
7, a sphere (there is, of course, a zero order in (k x ) 
term, but this obviously vanishes upon integration over 
the Fermi surface so we have omitted it). As the in- 
tegration is carried over a sphere, it is obvious that 
I(kf) 2 = f(k F ) 2 = (fc F ) 2 /3 and Jk F k F = 0. There- 
fore (B x ) = 0, and no spin splitting is expected from 
this procedure, even though in the Dresselhauss term 7 
maybe finite. Note that this cancellation would not hap- 
pen if the spin-orbit coupling term were /c-linear since 
the derivative of Bi nt would just be a constant. This is 
exactly what happens in the strained samples of GaAs 
where the spin splitting was explained by fc-linear terms 

A related fact shows that, due to its symmetry, the 
bulk-Dresselhaus term produces no uniform magnetiza- 
tion in the bulk of the sample. This is an easily falsifiable 
prediction of our theory. The Hamiltonian J5J) has two 
energy levels E± = |^fc 2 ± B where B = y/B~B~. The 
uniform magnetization (<jj) induced by an electric current 
Jj = dH/dkj (due to the applied electric field Ej) can 
be easily computed in linear response, and one obtains: 



Qij = (To-iJj 
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where np ± are the Fermi functions of the two bands. 
By inspection, all the components of Bi^- — are 
odd in the components fcj and hence vanish under inte- 
gration due to cubic symmetry. This leads to (a 1 ) = 0. 
By contrast, a fc-linear internal magnetic field, as in the 
strained samples, gives a finite uniform magnetization 
due to the fact that is a constant while B is isotropic 
in (and proportional to) k 01- The bulk Dresselhauss 
term is most likely also the explanation of the contradic- 
tion between the observed spin splitting along the [110] 
and [110] directions, and the uniform magnetization on 
these directions. In it is observed that although the 
spin splitting for S||[110] is consistently larger than the 
spin splitting for the £7| | [110] , the uniform magnetiza- 
tion for £:|| [110] is usually lower than that for £7||[110]. 
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This can be very well explained by the missing bulk- 
Dresselhauss term in the case when this term subtracts 
from the splitting on the [110] direction but adds to the 
splitting in the [110] directions, which makes it qualita- 
tively possible that the uniform magnetization observed 
be in agreement with the spin splitting. Quantitative 
modelling of this involves precise knowledge of the sam- 
ple Hamiltonian (including the strain) which is currently 
not possible (mainly because of limited knowledge about 
out-of-plane spin-orbit coupling terms). 

Even though it creates no uniform magnetization, the 
bulk Dresselhauss term does give rise to an intrinsic spin- 
Hall effect when under the action of an electric field. We 
define the spin current as usual (e = h 2 k 2 /2m): 



, 1 (dH A de dBi 



(9) 



and the expanded expression for the Green's function 
G(k,ioj n ) = [iuj n — H]" 1 as: 

G(k, iui n ) = f(k, iuj n )(g(k, iuj n ) + B^k)^) 
f(k, iuj n ) = ( ^_ £( i )) y_ B ^ g( k , «Wn) = - e(fcjjlO) 

When subjected to the action of an electric field E, the 
frequency dependent spin conductance (not including the 
vertex correction) can be found in linear response as: 



J 
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= vp Ek,n Tr[G(k, i(u> n + v m ))J\ (k)G(k, iu} n )Jj (HU) 

Summing over the Matsubara frequencies iu) n , analyti- 
cally continuing iv m — + to, as well as omitting a dissipa- 
tive term which vanishes upon momentum integration, 
we obtain the expression for the frequency dependent (re- 
active) spin conductivity: 
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where i,j,l,r,n = x,y,z. Unlike the uniform magneti- 
zation case, the integrand is even in k and finite upon 
integration. Hence the spin current in the unstrained 
GaAs can be qualitatively explained by the presence of 
a intrinsic spin-Hall effect due to the bulk-Dresselhauss 
term. Working in spherical coordinates, substituting the 
explicit expression for the bulk Dresselhauss spin split- 
ting Bi(k), using the expression (valid for small 7) of the 
difference between the Fermi momenta of the two spin- 
split bands: 



k F_,F^2mB(k F ). , F _ k F + k F 
- + ~ h 2 k F ' 2 ' 



(13) 



as well as integrating over the spherical angles, one ob- 
tains for the DC spin-Hall conductivity: 



k F 

\2ir 2tllJ 
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This is the intrinsic spin-Hall conductivity in the clean 
limit, H/r « B(k), which we call a 'intrinsic- For the car- 
rier concentration in pj we have k F = 10 8 to _1 and we 



obtain cf^^insic 



200fT 



" 1 . This is much larger than 



the observed conductivity of 0.2 ~ 0.5S1 



But this 



is expected since we have so far not taken the influence 
of disorder into account . In the experiment by Kato 
et al, the lifetime broadening due to impurity scatter- 
ing is much larger than the weak spin splitting due to 
the Dresselhauss coupling. The intrinsic spin-Hall effect 
is therefore in the dirty limit, and a significant reduction 
from the clean result is therefore expected. Note that the 
lifetime broadening due to impurity scattering leads to a 
reduction of the intrinsic spin-Hall conductivity. This is 
different from the extrinsic spin-Hall effect due to spin- 
orbit coupling to the impurity potential, which makes a 
small, but positive contribution to the spin-Hall conduc- 
tivity. 

To properly take into account disorder, one must per- 
form a self-consistent calculation taking into account 
both the self-energy and the vertex correction. In the 
case of the electron Rashba model with a fc-linear spin 
splitting, many groups have shown that the vertex cor- 
rection cancels the intrinsic spin- Hall effect BE! . How- 
ever, this cancellation seems to be special to the fc-linear 
spin splitting, and it has been shown that the vertex cor- 
rection due to fc 2 light/heavy hole splitting in the Lut- 
tinger model, or due to fc spin splitting in the heavy hole 
band, vanishes identically [8l|2lJ. We expect that for a 
similar reason, the vertex correction due to the fc 3 Dres- 
selhauss spin splitting would not cancel the intrinsic spin- 
Hall effect either. We will hence neglect the vertex correc- 
tion and focus on the self-energy correction which is easy 
to extract analytically. The self-energy approximation 
to disorder can be simulated by letting w = ih/r in Eq. 
(JT2J. The values in [lj are in the regime h/r >> B(k), 
we thus obtain: 

h 2 1 f d 3 k n E _ - n E+ dB r 

~ ^i^lnr^n 



2m{h/T) 2 J (2tt) 3 
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which, upon momentum integration gives the lower 
bound for the spin conductivity: 



7fc| 
105tt 2 \H/t 



(16) 
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For the values 7 = 25eVA , kp = 10 to we ob- 
tain a bulk Dresselhauss spin splitting energy 'ykp « 
O.Q25meV while h/r as 1.6meV for a sample of mobility 
fi = lm 2 /Vs as the one in the experiment. Using these 
values, we obtain for the intrinsic, disorder quenched 



spin conductivity a 



dirty 
intrinsic 



0.02O- 



This lower 



bound is smaller than the measured conductivity (which 



is 0.21T 



for small electric field and 0.5J7 



for 



large electric field. This is a lower bound for the spin con- 
ductivity since h/r is an upper bound for the frequency u> 
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in the dirty limit. Considering the uncertainty associated 
with the value of 7, the crudeness of the estimate, and 
the indirect determination of the experimental value, the 
agreement is reasonably good. 

The application of strain induces two extra spin split- 
tings in the Hamiltonian which are linear in the momen- 
tum k 17] . There is one structural inversion asymmetry 
(SIA) splitting of the form a(k y a x — k x a y ), and a bulk- 
inversion asymmetry (BIA) of the form b(k x a x — k y a y ), 
where a and 5 are strain dependent. For the values of the 
splitting in sample E used in we have a/H = 183m/s 
and 5/h = 112m/ s. We observe that the splitting at 
the Fermi momentum is 0. Oilmen for the SIA term and 
0.007meV for the BIA term. By contrast, the k 3 Dres- 
selhauss coupling is 0.025meV, so it is likely that it will 
dominate (although not overwhelmingly) the spin cur- 
rent. Moreover, a vertex correction computation for the 
SIA or BIA term separately reveals that the spin current 
caused by these terms vanishes upon the introduction of 
impurities |l9l | (exact numerical diagonalization results 
22, 2jj are, however, at odds with [ijj), whereas a vertex 
calculation for a fc 3 term shows finite spin current |sLl2*l|. 
It is therefore plausible that the bulk Dresselhauss term 
dominates the spin-Hall transport even in the strained 
samples. This naturally explains the independence of 
spin current on the crystallographic directions of the ap- 
plied electric field, since the bulk spin conductivity for the 
Dresselhauss term is direction independent. The exper- 
imental features observed can therefore be qualitatively 
explained by an intrinsic mechanism. 

In conclusion, we have shown that without any spin 
splitting in the electron band, the extrinsic spin-Hall ef- 
fect is far too small to explain the experimentally ob- 
served value of spin-Hall conductivity in In order to 
definitively determine the origin of the spin-Hall effect, 
we propose to carry out similar experiments in materials 
without any known intrinsic spin-orbit coupling, and a 
null result would give the definitive proof that the extrin- 
sic spin-Hall effect is far below the current experimental 
sensitivity, and can not be the origin of the spin-Hall 
effect observed in Rcf. [j. We have shown that the ex- 
perimental results are consistent with the interpretation 
of an intrinsic spin-Hall effect in terms of a bulk Dressel- 
hauss term in the unstrained sample. Furthermore, this 
intrinsic spin-orbit coupling is consistent with the appar- 
ent absence of spin splitting observedfin the unstrained 
samples) in the spin drag experiment We also predict 
that the uniform magnetization in the unstrained bulk 
GaAs samples will be close to zero due to the symme- 
try of the fc 3 Dresselhaus term. Since the experiment 



is carried out in a regime where the lifetime broaden- 
ing due to impurity scattering is large compared to the 
spin splitting, the observed spin-Hall conductivity is sig- 
nificantly reduced from the value in the clean limit. It 
is argued that the /c-linear terms in the strained GaAs 
samples 0, , although crucial for the appearance of a 
uniform magnetization, have a limited effect on the spin 
current due to the dominance of the Dresselhauss term, 
thereby qualitatively explaining the direction indepen- 
dence of the spin-Hall effect observed in the experiment. 
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